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We discuss a cosmology in which cold dark-matter particles decay into relativistic particles. We 
argue that such decays could lead naturally to a bulk viscosity in the cosmic fluid. For decay 
lifetimes comparable to the present bubble age, this bulk viscosity enters the cosmic energy equation 
as an effective negative pressure. We investigate whether this negative pressure is of sufficient 
magnitude to account fo the observed cosmic acceleration. We show that a single decaying species 
in a A = 0, flat, dark-matter dominated cosmology can not reproduce the observed magnitude- 
redshift relation from Type la supernovae. However, a delayed bulk viscosity, possibly due to a 
cascade of decaying particles may be able to account for a signiflcant fraction of the apparent cosmic 
acceleration. Possible candidate nonrelativistic particles for this scenario include sterile neutrinos 
or gauge-mediated decaying supersymmetric particles. 

PACS numbers: 98.80.Cq, 98.80.-k, 95.30.Cq, 95.35.-|-d, 95. 36. -fx, 98.80.Es, 98.62.Py 



I. INTRODUCTION 

A significant challenge facing modern cosmology is that 
of understanding the nature and origin of both the dark 
energy responsible for the present apparent acceleration 
and the dark matter j2| responsible for most of the 
gravitational mass of galaxies and clusters. The simplest 
particle physics explanations for the dark matter is, per- 
haps, that of a weakly interacting massive particle such 
as the lightest supersymmetric particle, an axion, or an 
electroweak singlet (e.g. "sterile" neutrino). The dark 
energy, on the other hand is generally attributed to a 
cosmological constant, a vacuum energy in the form of 
a "quintessence" scalar field possibly very slowly evolv- 
ing along an effective potential, or even relativistic effects 
derived from the deviation of the present matter distribu- 
tion from Friedmann homogeneity Q- See Q for recent 
review. 

In addition to these explanations, however, the sim- 
ple coincidence that both of these unknown entities cur- 
rently contribute comparable mass energy toward the clo- 
sure of the universe begs the question as to whether they 
could be different manifestations of the same physical 
phenomenon. Indeed, many suggestions along this line 
have been made for so-called unified dark-matter. One 
possibility is a dark matter composed of a generalized 
Chaplygin gas Q for which pressure depends upon den- 
sity p = —A/p"-. Although, it has been shown G] that a 
generalized Chaplygin gas produces an exponential blow 
up of the matter power spectrum which is inconsistent 
with observations, there are also more exotic proposals 
such as the flow of dark matter from a higher dimension 
0, or that the quintessence field itself can act as dark 



matter as in the Born-Infeld |8| model. 

The possibility of particular interest for the present 
work, however jis that of a bulk viscosity within the cos- 
mic fiuid (e.g. |2])- Such a term resists the cosmic expan- 
sion and therefore acts as a negative pressure. Indeed, 
it has been shown that for the right viscosity coeffi- 
cient, an accelerating cosmology can be achieved without 
the need for a cosmological constant. 

Although cosmic bulk viscosity is a viable candidate for 
dark-energy, to date there has been no suggestion of how 
it could originate from known physics and known particle 
properties. In this paper we consider a simple mechanism 
for the formation of bulk viscosity by the decay of a dark 
matter particle into relativistic products. Such decays 
heat the cosmic fluid and lead to an increase in entropy 
and are inherently dissipative in nature. Moreover, they 
lead to a cosmic fluid which is out of pressure and tem- 
perature equilibrium and can therefore be represented by 
a bulk viscosity. We propose a form for this viscosity and 
show that decay lifetimes comparable to the present bub- 
ble time naturally produce an accelerating cosmology in 
the present epoch. 

In the next section we summarize the general form for 
the bulk viscosity. Following that we consider its affect 
on cosmology and suggest a form for the bulk viscos- 
ity induced by particle decay. In the following section 
we discuss constraints on the properties of such parti- 
cles and argue that several candidates exist. In Section 
0we compute the magnitude-redshift relation for Type 
la supernovae in this cosmology and show that a single 
decay does not reproduce these data. Only if decays are 
delayed, e.g. by a cascade of particle decays, can a sig- 
nificant fraction of the cosmic acceleration be explained 
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in this scenario. 



The existence of such dissipation leads to a modification 
of the perfect-fluid energy-momentum tenor, 



II. BULK VISCOSITY FROM DECAYING 
DARK MATTER 

The fundamental problem that we want to address is 
the effect of the decay of non-relativistic dark-matter par- 
ticles into relativistic neutrinos. Bulk viscosity is a way 
to introduce the effects of this decay on the equations for 
cosmic expansion. The physical origin of bulk viscosity 
in a system can be traced to deviations from local ther- 
modynamic equilibrium. This can be illustrated with a 
simple abstract example. 

Suppose that the energy-momentum tensor in an ex- 
panding volume has contributions from both a compo- 
nent of particles obeying non-relativistic kinematics and 
a component following relativistic kinematics. Imagine 
that in a time step the system expands, but the momenta 
of the relativistic and non-relativistic particles redshift 
(i.e. change) differently. In effect, this causes these two 
components to have different "temperatures" describing 
their energy-momentum distribution functions. 

The second law of thermodynamics tells us 0| 
that the re-establishment of thermal equilibrium through 
(particle decay or) scattering of these component parti- 
cles off each other or on another medium is a dissipative 
process that will generate entropy. This entropy gener- 
ation can be related to the expansion rate or the local 
fluid velocity through a bulk viscosity term. 

Thus, bulk viscosity arises any time a fluid expands 
to rapidly and ceases to be in thermodynamic equilib- 
rium . The bulk viscosity, therefore, is a measure of the 
pressure required to restore equilibrium to a compressed 
or expanding system [l^ IT^ 120]. Hence, it is natural 
for such a term to exist in the cosmologically expanding 
universe anytime the fluid is out of equilibrium. Usually, 
in cosmology the restoration processes are taken to be 
so rapid that the establishment of equilibrium is almost 
immediate. However, there is a finite time for the system 
to adjust to the change of the equation of state induced 
by particle decays. For the cosmology proposed here, 
the attainment of equilibrium as the universe expands is 
delayed by the gradual decay of one or more species to 
another which occurs over ^ 10^" yrs. This leads to non- 
trivial dependence of pressure on density as the universe 
expands, and therefore a bulk viscosity. 

To see how this enters quantitatively in cosmology, we 
begin by summarizing the general treatment of imperfect 
fluids of Weinberg 15]. It will provide further insight into 
the nature of the bulk viscosity. 

When a fluid expands (or is compressed) and departs 
from thermodynamic equilibrium the processes that re- 
store equilibrium are irreversible. Hence, they are in 
general accompanied by an increase in entropy which is 
evidenced in the dissipation of energy. For the case of 
interest here, the increase in entropy and dissipation is 
the heating and pressure produced by the particle decays. 



T''" = ip + ppf'U, 



(1) 



where p and p denote density and pressure while 
is the four velocity. Processes of heat flow and shear 
can play no role in the Friedmann-Lemaitre-Robertson- 
Walker (FRLW) homogeneous and isotropic conditions 
of interest here. Hence, the only possible non-adiabatic 
dissipative contribution AT^'^ which guarantees transla- 
tional and rotational invariance for a fluid in motion with 
four velocity U" is given by ,15.] 
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(2) 



where a is the cosmic scale factor as specified below and 
C is the bulk viscosity coefflcient. The total energy- 
momentum tensor is 



(3) 



From Eq. ^ it is obvious that the effect of bulk viscosity 
is to replace the fluid pressure with an effective pressure 
given by. 



Peff =P-C3- 



(4) 



Thus, for large C it is possible for the negative pressure 
term to dominate and an accelerating cosmology to en- 
sue. It is necessary, therefore, to clearly deflne the bulk 
viscosity for the system of interest. 



III. COSMOLOGY WITH BULK VISCOSITY 

To examine the effect of the bulk viscosity from particle 
decay on the cosmic acceleration, we analyze a flat (k — 
0, A = 0) cosmology in a comoving FLRW metric. 



^2 gjj^2 fljj,2 



gf^^dx^'dx" = -dt^ + a^(t) 

for which U° ^ 1, [/' 0, and U)^ = 3a/a. 

We consider a fluid with total mass-energy density p 
given by. 



P = Pdm + Ph + Ph + Pj + Pi 



(6) 



where ph is the baryon density, pdm is the contribution 
from stable dark matter, ph is the density in unstable 
decaying dark matter, pi is the produced relativistic en- 
ergy density from decay while pj is any other relativistic 
matter, i.e. photons and neutrinos from the big bang. 
Because of decay, neither the total energy density in rel- 
ativistic particles pi = p^ + pi nor the pressure p = /9r/3 
is negligible for this cosmology even at the present epoch. 
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In the FLRW frame, the energy-momentum tensor 
(Eq. I2J) then reduces to 



Tqo 
T ■ 



P- 3C- ]9ij 



(7) 
(8) 

(9) 



where again, this last equation shows that the bulk vis- 
cosity enters as an effective negative pressure (i.e. dark 
energy) in the energy-momentum tensor. 

The Friedmann equation does not depend upon the 
effective pressure and is exactly the same as for a non- 
dissipative cosmology, i.e. 



-nGp 



(10) 



where p is the total mass-energy density from matter and 
relativistic particles (Eq. inj. 

Although absent from the Friedmann equation, the 
bulk viscosity does appear in the conservation condition 
Tt^'f^ = 0. To illustrate this consider a flat fc = 0, A = 
cosmology and ignore the small contribution from pb and 
initial background radiation p^. 



P = Ph+ Pi 



(11) 



The conservation equations can be solved to give the 
energy densities in matter and radiation: 



1 



Ph 



rPmoe 



-t/T 



and 



PI = -T 



Pia 



^ f e-''/^a{t')dt' + PBV 
^ Jo 



(12) 



(13) 



where pBV is the dissipated energy in light relativistic 
species due to the cosmic bulk viscosity, 



PBV ^9 C(i') - ait')Ut 



(14) 



The total density for the Friedmann equation will then 
include not only terms from heavy and light dark matter, 
but a dissipated energy density in bulk viscosity. This is 
the term that contributes to the cosmic acceleration. 



A. Bulk Viscosity Coefficient 

Bulk viscosity can be thought of 0,0,0] as a relax- 
ation phenomenon. It derives from the fact that the fluid 
requires time to restore its equilibrium pressure from a 
departure which occurs during expansion. The viscos- 
ity coefflcient ^ depends upon the difference between the 
pressure p of a fluid being compressed or expanded and 



the pressure p of a constant volume system in equilib- 
rium. Of the several formulations the basic non- 
equilibrium method is identical |15|| with Eq. Q . 



C3- = Ap 



(15) 



where Ap = p — p is the difference between the constant 
volume equilibrium pressure and the actual fluid pres- 
sure. 

In Ref. TH'l the bulk viscosity coefficient is derived for a 
gas in thermodynamic equilibrium at a temperature Tm 
into which radiation is injected with a temperature T and 
a mean thermal equilibration time t^. The solution for 
the relativistic transport equation [l8j can then be used 
to deduce infer [l^ the bulk viscosity coefficient. For 
this case the form of the pressure deficit and associated 
bulk viscosity can be deduced from Eq. (2.31) of Ref. 0| 
which we modify slightly and write as. 



Ap . 



T) 



Ap. 
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(16) 

where the factor of 4 comes from the derivative of the 
radiation pressure p ~ of the injected gas, and the 
term in brackets derives from the detailed solution to 
the linearized relativistic transport equation [jj. This 
term guarantees that no bulk viscosity can exist for a 
completely relativistic gas. In the cosmic fluid, however, 
we must consider a total mass-energy density p given by 
both nonrelativistic and relativistic components. 

In the present context we also have a thermalized gas 
into which relativistic particles at some effective temper- 
ature are being injected. The deflcit from equilibrium 
pressure, however, is due to the presence of unstable de- 
caying nonrelativistic dark matter. At any time in the 
cosmic expansion the pressure deficit will be 1/3 of the 
remaining mass-energy density of unstable heavy parti- 
cles. Hence, we replace p^/3 with ph/3 in Eq. and 
write, 



Ap = 



4phTo 



5^ 



(17) 



Here, the equilibration time Te is determined (16j from 
the particle decay time r. 



Apjt) 
Ap{0) 



dt = 



[I - 3(d/a)T] 



(18) 



where Ap(0) denotes the initial pressure and the denom- 
inator results from approximating H — a/ a constant. 
Note, that this factor acts as a limiter to prevent unreal- 
istically large bulk viscosity in the limit of a large r. 

Following the derivation in J5] , and inserting Eq. I|17() 
in place of Eq. 116|l . we infer the following ansatz for the 
bulk viscosity of the cosmic fluid due to particle decay, 



1 - 



P\ + Pi 



(19) 
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where the square of the term in brackets comes from 
inserting Eq. H17II into the hnearized relativistic trans- 
port equation of Ref. |0| . Equation (|19|l implies a non- 
vanishing bulk viscosity even in the limit of long times 
as long as the total mass energy density is comprised 
of a mixture of relativistic and nonrelativistic particles. 
Hence, one should be cautious about using this linearized 
approximation in the long lifetime limit. Even so, a more 
general derivation has been made ^3 which shows that, 
even in the limit of interest here of a long radiation equi- 
libration time there is a non- vanishing bulk viscosity con- 
sistent with experimental determinations. 



IV. DECAYING DARK MATTER CANDIDATES 

Having postulated the existence of a decaying dark 
matter particle, it is important to briefly examine the 
constraints on such decays and whether such candidate 
particles could exist. To avoid observational constraints 
the decay products must have very little energy in pho- 
tons or charged particles. The implied background in 
energetic photons with an energy density comparable to 
the present matter energy density would have been easily 
detectable. Hence, the decay products must be in some 
form which is not easily detectable. Neutrinos would be 
a natural candidate for such a background. In this case 
there are several decaying dark matter possibilities which 
come to mind. 



A. Sterile Neutrinos 

Perhaps the most realistic possibility is the decay of a 
sterile neutrino into light "active" neutrinos [l^l- Mod- 
els have been proposed in which singlet "sterile" neutri- 
nos Vs which mix in vacuum with active neutrinos (uc, 
Vu.^ Vt\ provide warm and cold dark matter candidates 
[ialillilEllllll^ii- In most of these models the 
sterile neutrinos are produced in the very early universe 
through active neutrino scattering-induced de-coherence. 
This process could be augmented by medium enhance- 
ment stemming from a significant lepton number. In 
these sterile neutrino production processes there are two 
principal parameters: (1) the sterile neutrino mass mg; 
and (2) the sterile neutrino's vacuum mixing angle 9 with 
one or more of the active neutrino flavors. The net lep- 
ton number(s) of the universe could be regarded as an 
additional parameter. By virtue of the mixing with ac- 
tive neutrino species, the sterile neutrinos are not truly 
"sterile" and, as a result, can decay. For ra^ < lOMeV 
the dominant decay mode is into light, active neutrino 
species. The rate for this process is ,22j| 

. 91 IS /sin^2^?\/ ms \5 , , 

r.«(8.7xl0-s-)(^)(^) . (20) 



Likewise, there is a sub-dominant Vs decay branch into a 
light active neutrino and a photon with rate 

„ / oi T\ f sin^20\ / \^ , , 

r^.. « (6.8 X 10^23 -i\ I ^ \ 

^ ^ \ 10-15 y VlMeV/ ^ ' 

In this process the photon will be mono-energetic with 
an energy which is half the Vg rest mass. Because the 
primary i^g decay mode and the radiative branch scale 
he same way with rus and sin^ 29, there is a fixed ratio 
of these rates. The best particle candidates for a decay- 
induced bulk viscosity are those with a lifetime of order 
the Hubble time Hq^ and rest masses ~ 1 MeV. Setting 
= Hq we find that the relation between the Vs rest 
mass and vacuum mixing angle is 

,„..3.1Mev(— ) (--^) , (22) 

where h is the Hubble parameter at the current epoch in 
units of 100kms~i Mpc~^ and we have scaled our result 
toh^ 0.71, the WMAP best fit value. We conclude that 
one or a number of sterile neutrinos with rest masses 
in the ~ MeV range could provide a significant decay- 
induced bulk viscosity. 

Regarding observational constraints, let us note that 
our bulk viscosity-selected range for mg from Eq. H22|) is 
relatively insensitive to the vacuum mixing angle. How- 
ever, the radiative decay branch rate T^-y is linearly pro- 
portional to sin^ 29. Keeping mg 1 MeV, we can ad- 
just sin^ 29 so that the diffuse deca y p hoton flux is just 
at or below the observational limit '2^ |2^ Isof from the 
Diffuse Extragalactic Background Radiation (DEBRA). 
For this mg = 1 MeV the DEBRA limit would require 
sin^ 29 < 10-1^ 

We conclude that it is possible to meet the bulk vis- 
cosity lifetime requirement and (barely) get under the 
DEBRA limit with sterile neutrinos as the decaying dark 
matter. We also note that sterile neutrinos with these 
parameters (tos ~ IMeV, sin^ 29 ~ lO"!^) could be pro- 
duced in the early universe in the requisite relic densities 
{i.e., near closure) only in scenarios with large le pton 
number (s) and medium-enhanced de-coherence |22. l3l|. 
or with new neutrino couplings [2^ . 



B. Decaying Supersymmetric Dark Matter 

For supersymmetric dark matter candidates. It is gen- 
erally assumed Q that the initially produced dark mat- 
ter relic must be a superWIMP in order to produce the 
correct relic density. Later, this superWIMP is then pre- 
sumed to decay to a lighter stable dark matter particle. 
One interpretation of such a candidate for decaying dark 
matter, is then a decaying superwimp with a lifetime 
comparable to the present bubble time. 

Alternatively, the light supersymmetric particle itself, 
might a candidate for decay. If the dark matter is a light 
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unstable supersymmetric particle, then one might imag- 
ine an R-parity violating decay. In one scenario a par- 
ticle might decay by coupling to right-handed neutrinos 
which then decay to normal neutrinos. Another possi- 
bility could be gauge-mediated supersymmetry breaking 
involving the decay of a supersymmetric sneutrino into a 
gravitino plus a light neutrino. 



V. RESULTS 

Having defined the cosmology of interest we now ex- 
amine the magnitude-redshift relation for Type la super- 
novae (SNIa). The apparent brightness of the Type la 
supernova standard candle with redshift is given |32| by 
a simple relation for a flat A = cosmology. The lumi- 
nosity distance becomes, 



-1/2 



,(23) 



where Hq is the present Hubble parameter. The Qi are 
the energy densities normalized by the critical density at 
each epoch, i.e. ^li{z) = ^ttG pi{z) / iH^ . ilh is the clo- 
sure contribution of the decaying heavy cold dark matter 
particle which is taken to produce light neutrinos fii or 
other relativistic particles as it decays. Note that fih, 
and U,\ each have a nontrivial redshift dependence due 
to particle decays, while stable dark matter and baryons 
^'dm{z) + ^h{z) obey the usual (1 -f z)~'^ dependence 
with redshift. Here, and in the following discussion we 
will define VIm as the present sum of nonrelativistic mat- 
ter, i.e. JIm = f^h(-z = 0) + ilDM(z = 0) -I- VL\,{z = 0.) 



(Gyr) 



0,M 0.A Xr 



00 (ACDM) 0.31 0.69 1.14 
20 0.16 0. 3.89 

1 X 6, 20 0.37 0. 1.90 
oo (CDM) 1. 0. 3.23 



TABLE I: Parameter sets for various fits to the SNa 
luminosity-redshift relation for Hq = 71 km s~^ Mpc~^ and 
rib = 0.044. In the decaying (finite r) models no stable dark 
matter was assumed (i.e. f^DM = 0). 

Figure ^ compares various cosmological models with 
some of the recent combined data from the High-Z Super- 
nova Search Team and the Supernova Cosmology Project 
0, l33l | , while Table 1 summarizes the relevant parame- 
ters and reduced goodness of fit. The lower figure 
shows the K-corrected magnitudes m = 7\/ + 5 log£'i + 25 
vs. redshift plotted relative to an open flmiDAi, ^b, = 
0, fife = 1 cosmology. 

The solid line on the upper and lower graphs in Figure 
n shows the result of adding bulk viscosity from particle 



decay. The upper figure gives the distance-redshift rela- 
tion while the lower figure shows the evolution of mag- 
nitudes relative to a fiducial fifc — l/(aoi?o)^ — 1 open 
cosmology, for which 



DLi^k = 1) 



2Ho 



z + 1- 



1 



(^ + 1) 



(24) 



and the relative distance modulus is given in the usual 
way A(m - M) = 5 log = 1)]. 

From the lower graph of Figure^we see that, although 
the bulk viscosity has indeed provided a negative pressure 
it does not reproduce the supernova distance-red shift re- 
lation. In fact it is much worse than the usual ACDM 
cosmology and is even worse than a pure matter domi- 
nated cosmology. The reason for this can be discerned 
from Figure |21 Although the bulk viscosity is substan- 
tial, it scales with the decaying dark matter which falls 
off faster with time than because of the decay. An 
accelerating cosmology requires a nearly constant value 
of ptot with time. 

A flattening of ptot could be achieved in this context 
if the onset of the bulk viscosity could be delayed until 
near the present epoch due to a cascading decay. In this 
possibility, the flnal decay would produce the relativistic 
products and the bulk viscosity. This final decay, however 
is preceded by a series of decays to nearly degenerate 
states with shorter lifetime. 

The cascade possibility might occur, for example, 
among sterile neutrinos. Another cascade possibility Q 
is that the that the initially produced dark matter relic is 
a super WIMP. Later this super WIMP decays to a lighter 
stable dark matter particle, e.g. a gravitino. One could 
have a cascade of decaying superWIMP states to a fi- 
nal unstable state or to unstable light supersymmetric 
particle states. 

Figure O illustrates a possible evolution of energy den- 
sity in this scenario. Here, a decay among six states each 
with a lifetime of ti 1 Gyr each is followed by the decay 
of a long lived particle with r = 20 Gyr for the case in 
which all matter starts in the first member of the cascade. 
The activity in the final decay product is delayed by the 
time needed to decay through the intervening states. The 
rates of the final and intermediate decays decay are given 
by a solution to the Bateman equation whereby the abun- 
dance of the final product is given by 



Ph — ^hj exp {—Xjt) 



(25) 



where = Tj ^ is the decay rate of each species and the 
hj are given by. 



(A. -A,) 



(26) 



For this possibility, we found that it was not possible 
in this way to completely account for the cosmic acceler- 
ation, though a significant fraction could be obtained as 
illustrated by the dash-dot-dot lines on Figure ^ 
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VI. CONCLUSION 

We have considered models in which the apparent cos- 
mic acceleration is affected by the bulk viscosity pro- 
duced from the decay of a dark matter particle to light 
relativistic species. An expression for the bulk viscosity 
is deduced and the implied redshift-distance relation has 
been computed. 

As an illustrative example we considered the decay of 
dark matter with a lifetime of 20 Gyr in this cosmology. 
From the reduced Xr values in Table 1, and the lines in 
Figure H it is apparent that a flat A = cosmology with 
bulk viscosity from decay of a single dark-matter species 
does not do better than a ACDM or a matter dominated 
cosmology. This is because the total mass-energy den- 
sity does not become nearly constant with scale factor, 
but falls off more rapidly than even a simple matter domi- 
nated cosmology due to the combined effects of the decay 
of the dark matter and the emergence of a high density 
of relativistic particles. We show, however that if the 
emergence of the bulk viscosity is delayed, then some, 
but not all, of the acceleration required by observations 
of type la supernovae at high redshift can be explained. 
As we have outlined above, one mechanism for delaying 
the bulk viscosity could be a cascading decay process. 

Obviously, however, one must decide whether the 
dilemma of a cosmological constant is less plausible than 
the dilemma of bulk viscosity produced by a delayed cas- 
cade of decaying dark matter particles. Our goal here, 
however, has merely been to argue that the possibility ex- 
ists. Having established that at least a possible paradigm 
exists, in future work we will examine the possible influ- 
ence of this scenario on the CMB and the growth of large 
scale structure which will also constrain this possibility. 



Indeed, a number of recent studies (e.g. 0) suggest 
that changes in the extrinsic curvature due to changing 
relativistic gravity in an inhomogeneous cosmology can 
lead to cosmic acceleration. We are currently developing 
a computer model similar to |34| that will include the 
decay of heavy neutrinos. Our preliminary analysis indi- 
cates that the cosmic acceleration in a universe with this 
type of dark matter decay may be enhanced by its effect 
on large scale structure, i.e. decay of heavy neutrinos in 
a non homogeneous cosmos may increase the expansion 
effect. 

In brief, the decay produces a flow of light neutrinos 
from galactic clusters. Given a decay time r and galactic 
clusters separated by a distance L this flow produces a 
momentum density of the order s = pL/r (in units of 
c — G = V). From the momentum constraint [Eq. (A4) 
of 34]], an enhancement of the extrinsic curvature of or- 
der pL^/r will occur. From the Hamiltonian constraint 
[Eq. (A3) of [S^]] the trace of the extrinsic curvature will 
be reduced by a factor the order of {6K/K)^ and since 
a/a ^ K /Z, a/ a becomes more nearly constant implying 
acceleration. 
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FIG. 1: Evolution of luminosity distance with redshift for a 
cosmology with bulk viscosity. Points are from the Gold data 
set of 33]. The upper figure shows the luminosity distance 
vs. redshift. The lower figure shows the evolution of magni- 
tudes relative to a fiducial Qk ~ 1 open cosmology. In each 
figure the upper dashed line shows the evolution of a stan- 
dard ACDM cosmology and lower dot-dashed line shows the 
evolution of an = 1 cosmology. The solid line is for a 
illustrative decaying dark matter model with t — 20 Gyr. the 
dash-dot-dot line illustrates the evolution of a cosmology in 
which a cascade of six particle decays each with a lifetime of 
ri = 1 Gyr is followed by a final radiative decay with r = 20 
Gyr. 
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FIG. 2; Evolution of the quantities p-y, pM, Pbv, and ptot as 
labeled for a cosmology in which the dark matter decays with 
a lifetime of 20 Gyr 
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FIG. 3: Evolution of the quantities p^, pM, pBV, and ptot 
as labeled in a cosmology in which decays among six nearly 
degenerate states occur with a lifetime of t\ 1 Gyr each is 
followed by the decay of a long lived particle with r = 20 
Gyr. 



